GAUSSIAN RANDOM MATRIX MODELS FOR q DEFORMED 
GAUSSIAN VARIABLES 



PIOTR SNIADY 



Abstract. We construct a family of random matrix models for the q— 
deformed Gaussian random variables G M = a M + a* where the annihilation 
operators a M and creation operators a* fulfil the g-deformed commutation re- 
lation a^a* — qa^a^ = r M „, r is the covariance and < q < 1 is a given 
number. Important feature of considered random matrices is that the joint 
distribution of their entries is Gaussian. 



1. Introduction 

1.1. The deformed Gaussian variables. The g-deformed Gaussian random 
variables G M — + a*, where operators a M and their adjoints a* fulfil deformed 
commutation relations 



(1) 



were introduced by Bourret and Frisch | FB | . These operators act on a Hilbert space 
K. which has a unital vector f2, called a vacuum, with the property that 



(2) 



for every value of the index fi. 

With the help of the vector f2 one can introduce a state r on the algebra of 
operators acting on K. as follows t(X) = (ft, XCl). The state r plays a role of the 
non commutative expectation value. From (0) and (H) follows [ BS1 ] that for any 
m G N and any indexes . . . , [iim we have that 



(3) 



t(G Mi • • • G t _ l2m _ 1 ) — 0, 



(4) r(G Ml • • ■ G M2m ) = ]T q tM T Cldl ■ ■ ■ r Cmdm , 

IT 

where the sum is taken over all pair partitions ir — {{ci, o?i}, . . . , {c m , c? m }} of the 
set {1, . . . , 2m} and i(7r) is the number of crossings of the partition ir. For the 
reader's convenience we shall recall definitions of a pair partition and of its number 
of crossings in Sect. |[ 

From the quantum probability point of view all the information about non com- 
mutative random variables G M is encoded in their moments t(G Mi •••G flm ) and 
therefore Eq. (|) and (|) can be treated as an alternative definition of g-deformcd 
Gaussian variables G M . 
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1.1.1. Applications of deformed Gaussian variables. Eq. (^|) and ^ show that for 
q = 1 operators G M have the same moments as classical Gaussian variables with 
the mean zero and the covariance r^,,, what should explain why do we call G a 
deformed Gaussian variables. Eq. (|]) is for q — 1 called the canonical (or bosonic) 
commutation relation. 

For other special choices of the deformation parameter q variables G^ also have 



natural probabilistic interpretations [FB|, namely as increments of a dichotomic 
Markov process (for q = — 1) or as Wigner's large random matrices (for q = 0). 
Voiculescu [VI] has made a remarkable observation that for q = random variables 



G a are free semicircular elements (an analogue of independent Gaussian variables 



in the free probability theory of Voiculescu |V3|, VDN|). Eq. (m) is for q = — 1 







called the canonical anticommutation relation (or fermionic relation) and for q 
is called the free relation. 

Therefore it was natural to expect that the relations (||) and (Q) which are a 
simple generalisation of the three mentioned above: bosonic, fermionic and free 
cases would give rise to interesting probabilistic obj ects. 

Indeed, it was observed by Bozejko and Speicher | BS1 | that a related to Eq. ([!]) 
Brownian motion is a one component of an n-dimensional Brownian motion which 
is invariant under the quantum group S ^qU{n) of Woronowicz for < q < 1. 

Another application of g-deformed Gaussian variables, this time as generalised 
quantum statistics was proposed by Greenberg |Gr[ and Speicher [{5p2f. 



The existence of operators and a* fulfilling deformed commutation relations 



(0) was proven by Bozejk o and Speicher | BS2 |. Later it was proven by Bozejko 
Kummerer, and Speicher |BKS| that the von Neumann algebra generated by q- 
deformed Gaussian variables Gi,G2,... (—1 < q < 1) is a 1I\ factor. There 
are today many open questions concerning these factors, particularly if they are 
different from the free group factors. 

In this paper we present a natural probabilistic representation of the g-deformed 
Gaussian variables for all q € [0, 1] a s so me random matrices, what was one of the 
open questions posed in the paper [FB|. A remarkable property of our model is 
that the joint distribution of entries of our matrices is Gaussian. 

Recently a related problem of finding a random matrix model for so called q— 
deformed circular system was solved by Mingo and Nica [MN|. 



1.1.2. The covariance T nv . Indexes /i, v are elements of a certain set 9Jt. A neces- 
sary and sufficient condition for operators G a to exist is that the function T av is 



positive definite [ BS2j |, i.e 



for all ai,..., a n £ K and fix, . 
covariance functions are: 



. . , fi n € 9Jt. Typical examples of sets 9Jt and 



371 = N and T^j = 5 l j. For q — 1 we have that Gi,G 2 , 



. . is a sequence 
we have that 



of independent, standard Gaussian variables while for q 
Gi, G2, . • • is a sequence of free semicircular elements VDN]. 
371 = R+ and r t s = min(f, s). For q = 1 we have that Gt is a Brownian 
motion, for q — we have that Gt is a noncommutative stochastic process 
with free increments 
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• 9Jt is a real Hilbert space and the covariance is denned by the scalar product: 
T^i/i = (4>i ip)- The case OJl = £ 2 (R + ) is often used in the white noise calculus. 

1.1.3. The distribution of a deformed Gaussian variable. A distribution of a ran- 
dom variable corresponding to the bounded selfadjoint operator G is a measure v 
supported on the real line R such that 



r(G n ) = / x n dv{x) 

for all n G N. 

It can be shown | ]gz| that the distribution v q of a g-deformed Gaussian variable 
with the variance equal to 1 is given by a measure v q supported on the interval 

with a density 



oo 

v q {dx) = -^/T~q sinfl TT(1 - <f )|1 - q n e 2ie \ 2 



dx. 



where 



with 9 G [0, 7r] 



2 

: COSW 



1.1.4. Canonical commutation relations and Ito's formula. It is not merely an ac- 
cident that for q = 1 there is a correspondence between the commutation relation 
(P and Gaussian random variables. 

If we consider a probability space generated by a Brownian motion B(t) then 
every real random variable X with a finite second moment can be uniquely expressed 
as a series of iterated Ito integrals 

X = X^+jr f X^(t 1 ,... ) t i )dB(t 1 )---dB(t i ) ) 

i=l J0<tx<---<ti<oo 

where X^ G K and for each i G N we have that X^ : (R+) 1 — > K is a symmetric 
function of its i arguments. By the bosonic Fock space we call the space of sequences 
(l''')i>o such that : (R + ) 1 — > R is a symmetric function. 

The Fock space carries a structure of a Hilbert space with a scalar product 

(X,Y) = E[XY] ^x (0) y (0) + 



y 1 

^ n! 

n>l 



X W (*!,..., tn)^ M (h,...,t n )dh ■■■dt n . 



The Fock space representative of the random variable constantly equal to 1 will 
be denoted by Q. We have Q (0) = 1, fi (n) (*i, . . . ,t„) = for every n > 1. Note 
that for every random variable we have 

E[X] =E[X1] = 

It is often convenient to work with the Fock space representations than with 
random variables themselves. An interesting question is to determine the Fock 
space representation of a product [J Q 4>(t)dB(t)]X if the Fock space representation 
of X is given and the integral is taken in the Ito sense. It is a simple corollary from 
the Ito's theorem that the multiplication by / °° 4>(t)dB(t) can be expressed as a 
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sum of two operators acting on the Fock space: so called annihilation operator 
and its adjoint, creation operator a*r. 



<j}{t)dB{t)= ait> + al. 

o 

They have the following properties: 

(5) a^ay - a^a^ = (0,-0), 

(6) a^fl = 0. 

We see that the commutation relations (||) are equivalent to the statement of 
the Ito's theorem. Conversely, postulating commutation relations ([!]) is equivalent 
to saying that the non commutative stochastic process G t fulfils some non com- 
mutative Ito's formula. Such a stochastic calculus for ^-deformed operators was 



considered by the author |Snl] 



1.2. Random matrices. For a general introduction to random matrices and their 
applications in mathematics and physics we refer to the monographs of Mehta M , 



Hiai and Petz [HP| and to overview articles [Br, GMGW|. 

There are essentially two kinds of questions concerning eigenvalues of a random 
matrix one can ask. Global questions are of the type: what is the asymptotic dis- 
tribution of eigenvalues if the size of a matrix is large enough, while local questions 
concern, for example, the distribution of the spacings between consecutive eigenval- 
ues. The global questions are much easier to answer and free probability theory has 
provided powerful tools to answer such questions for many random matrix model 



1 [V| |VDN| , |b| 



Recently random matrices were used as a powerful tool in the theory of C*- 



algebras by Haagerup and Thorbjoersen |HT | 



1.3. Overview of the paper. In Sect. || (which is independent of the rest of this 
article) we present heuristic motivations for some matrix models considered in this 
article. 

In Sect. H we introduce the notations and construct a family we construct an 
auxiliary family of Gaussian random matrices 

In Sect. |we construct a central object of this paper, namely a family of Gaussian 
random matrices S^ 41 . Since the structure of matrices S is a bit complicated, it 
is convenient to think about them as some weighted sums of the auxiliary matrices 
R, which have a much easier structure. 

As the index N tends to infinity, the size of our matrices grows exponentially. 
We show that matrices S asymptotically have the same expectation values as q— 
deformed Gaussian random variables. 

Our construction bases on the observation that for a finite dimensional Hilbert 
space Ti. there are 2^ decompositions of a tensor power into two factors 7-t® N = 
n ®k u ®{N-k) which corr espond to ways of decomposing a set {1, . . . , N} = 
A U ({1, . . . , N} \ A) into two subsets. The appropriate isomorphisms ja ■ TL® N — ► 
H®\ A \ (g)H' 8 ( 7V -l A l) give rise to isomorphisms of matrix algebras ]a ■ End(H®' A ') <g> 
EndOtt^-l- 4 !)) -> End(H® N ). 

Each auxiliary matrix R( N ^' A '^ is obtained by embedding a small standard her- 
mitian random matrix R^ 41 ^ ^ End(7i'^') into a bigger algebra End(H Ar ). The 
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embedding is implemented by the isomorphism ja- We have that 

AQ{1,...,N} A*Z{l,...,N} 

where <j a N ^ is a certain weight function. 

It turns out that the the commutation properties of two random matrices 
ft(N),Ai,n corresponding to two decompositions given by sets Ai,A 2 depend on 
the number of common elements of Ai and A 2 . For example, if Ai n A 2 — then 
these two matrices commute and, informally speaking, the more elements A\ and 
A 2 have in common, the more they behave like a pair of free random variables. 
Therefore the expectation value of a product of many matrices R can be evaluated 
from the number of elements of Ai n Aj and in this way we are able to calculate 
the mixed moments of matrices S. 

The choice of a normed weight a A N ^ is equivalent to a choice of a probabilistic 
measure p on the set of all subsets of {1, . . . , N} defined such that the measure 
of a singleton {A} is equal to (o^f'*) 2 for any A C {l,...,iV}. If this measure 
is, loosely speaking, concentrated on the sets of order cs/~N then if A and B are 
independent random variables with distribution given by the measure p then |AnS| 
is asymptotically Poisson distributed with the parameter A = c 2 . For appropriate 
choice of c we are able in this way to obtain a random matrix model for g-deformed 
Gaussian random variables for all < q < 1. 

In Sect. |^ we show that our matrices converge to g-deformed Gaussian variables 
not only in the sense of expectation values mixed moments, but that that mix 
moments converge almost surely. 

Sect. ^ is devoted to proofs of technical lemmas. 

The random matrices S considered in this article are complex hermitian. How- 
ever there are no difficulties to extend these results to real symmetric or symplectic 
hermitian. 

This paper shows that the g-deformed probability theory, which was regarded 
until today as purely abstract and algebraic, in fact has natural probabilistic models 
just like the free probability theory has. 



2. Heuristics of Random Commutation Relations and Random 

Gaussian Matrices 

This section an an independent part of the article and notations used here will 
not be used in the subsequent considerations. We have to warn the reader that this 
section is very informal, however by such informal considerations it is much easier 
to get an insight to the nature of the problem. 



2.1. q— deformed Gaussian variables. Our motivations how to find random ma- 
trices which asymptotically behave like g-deformed Gaussian variables were inspired 
by a careful study of the article of Speicher [ 3pl | . 

In this paper he shows a certain non commutative central limit theorem that 
if a suitably normalised family of centered non commutative random variables 
K\ , K 2 , . . . has a property that each pair of them either commutes (with prob- 
ability p) or anticommutes (with probability 1 — p) and if for each pair the choice 
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of one these possibilities is made independently, then the distribution of a normal- 
ized mean Kl+ '^ Kn converges (as n — ► oo) to the distribution of a g-deformed 
Gaussian random variable with q = 2p — 1. 

Now we shall construct a family of random (non Gaussian) matrices which almost 
fulfils the assumptions of the Speicher's theorem. 

Let us fix a real number < q < 1. For any N <G N let us consider a family of 
2 N x 2 N matrices 

indexed by s E N, where for each s e N and 1 < n < N we have that K™ is a 2 x 2 
matrix chosen randomly according to the following table: 



matrix 


oo 


-co 


CTl 


-<Tl 


02 


-0"2 


0-3 


-0- 3 


probability 


l-'ir 
2 


l-3r 
2 


r 
2 


r 
2 


r 
2 


r 
2 


r 
2 


r 
2 



where < r < I is a real number to be specified later and 



oo 



" 1 


" 




' 


1 " 




' -i~ 




' 1 





1 


, o-i = 


1 





, 0"2 = 


i 


, o 3 = 


-1 



are Pauli matrices. The random choices of matrices K" should be made indepen- 
dently. 

These eight hcrmitian matrices ±<7j have a property that each pair of them either 
commutes or anticommutes. It is a simple calculation that the probability that two 
independent 2x2 matrices (chosen according to the above table) anticommute is 
equal to 6r 2 . 

Therefore all matrices K s are hermitian and furthermore each pair K s and K t 
(s 7^ t) either commutes (if the number of indexes n such that K™ anticommutes 
with K™ is even) or anticommutes (if this number is odd). We see that the proba- 
bility of the event that K s commutes with K t is equal to 

]T (^](l-6r 2 ) N - 2m (6r 2 ) 2m 



P 



while the probability that K s anticommutes with K t is equal to 

l-p= Y, ( 2 m 1 )(l-6r 2 ) JV -( 2m+1 )(6r 2 ) 2ro+1 . 

The difference of these two probabilities is equal to 2p — 1 = (1 — l2r 2 ) N , 

We would like to apply now the Speicher's theorem to the family (K s ) by choosing 
r as a function of N such that q = [1 — \2r 2 N ] N . For large TV one can take the 

approximate value — \J— jf^- 

Unfortunately the Speicher's theorem cannot be applied directly since it turns 
out that the events "K s commutes with K t " are not independent for different pairs 
{s,t}. However, if N tends to infinity it can be justified that they are, loosely 
speaking, more and more independent. 

By classical central limit theorem we have that the limit distribution of K = 
-j^{K\ + • • • + K n ) (regarded as a classical random variable with values in a vector 
space of matrices) is Gaussian. In order to characterise this distribution uniquely 
we have to give the mean and the covariance of entries. The mean value of each 
entry of the matrix K is equal to and the covariance factorises as follows 

E [- K '(ii,...,iiv),(ji,...,jjv)- K '(fei,.-->fejv)>(«i,--->«Jv)] = 
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- E [ K (il,...,iN),(h,-dN) K (h,-M,{k u ...,k N )] - 

E\K} t Kl , ]->-fflKF , K? , }. 

I fci ,t i J L in, 3n k N ,t N l 

Above we have used the fact that due to the factorisation M 2 n = Mi <g> • • • <£> M2 we 
can parametrise the rows and columns of a 2 N x2 N matrix by sequences (it, . . . , in), 
where < it, . . . , ijv < 1- 

A simple calculations shows that the covariance of 2 x 2 matrices given by the 
table above is equal to 

It turns out that this is exactly the covariance we will considered in Proposition 
[j] for the special case of d = 2, see (|l5|). 

3. Random Matrix Model for Deformed Gaussian Variables 

3.1. Pair partitions. 

Definition 1. A pair partition of a given finite set M is any decomposition of M 
into a family ir — {{ci, rfi}, . . . , {c m , d m }} of disjoint sets each having exactly two 
elements: 

Ci <k, for 1 < i < m, 
{ci,di}r\{cj,dj} =0, fori^j, 
{ci, di} U • • • U {c m , d m } = M. 

The sets {ci, di}, . . . , {c m , d m } are called lines of the pair partition ir. 

If M is an ordered set we say that two distinct lines {a, b},{u,v}, a < b, u < v 
cross ifa<u<b<v or u < a < v < b. 

For a given pair partition it we will denote by i(ir) the number of crossings in 
ir, i.e. number of all unordered pairs of lines {a,b},{u, v} £ ir such that the lines 
{a, 6}, {it, v} cross. 

Example. There is only one pair partition of a two element set {1, 2} and there 
are three pair partitions of a four element set {1,2,3,4}, namely {{1, 2}, {3, 4}}, 
{{1,4}, {2,3}}, {{1,3}, {2, 4}}. 

We have z ({{1,2}, {3, 4}}) = 0, i ({{1,3}, {2,4}}) =1. 

Sets having an odd number of elements do not have any pair partitions at all. 

3.2. Notations. Let us fix a natural number d > 2. For any N £ N we define 
H r = C d for 1 < r < N and H {N) = <S> 1<r<N U r . In the following we shall often 
omit the index (N) standing at various objects, however we have to remember 
about their dependence on N. 

If /o, . . . , fd-i is an orthonormal basis of C d then e- Y — fa ® • • • ® fi N is an 
orthonormal basis of 7i^ N \ where i = (it, ■ ■ ■ , «jv) and < it, ■ ■ ■ , iff < d — 1. 

Here and in the following by bold letters i, j, k, . . . we shall denote variables which 
index the basis of and always we have i = (it,..., In), k = (kt, ■ ■ ■ , fc/v), e ^ c - 

For any set A where A — {at, ■ ■ ■ , at] C {1, . . . , N}, at < ■ ■ ■ < a k and A' = 
{1,2,..., N} \ A = {bt, ■ ■ ■ , bN-k}, b\ < ■ ■ ■ < bN-k we consider Hilbert spaces 
K A = ® reA H r and H A ' = ® reA , H r and an isomorphism j A : H (N) ^ H A ® H A ' 
given by a grouping of factors: 

Vt ® •• • ® V N 1-^ (v ai ® • •• ® »o fc ) ® (Vb! ® • • • ® v bN _ k ). 
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This isomorphism induces an isomorphism of matrices ]a '■ End(7i' 4 )(E)End(7i" 4 ) — > 
End(?<W). 

In the following we shall denote by tr the normalised trace on End(7i ( - Ar - ) ) defined 
by tr(M) = ^Tr(M), where Tr denotes the standard trace. 

3.3. Iverson's notation. In the following we shall use sometimes Iverson's nota- 



tion [GKP] as an alternative to the Kronecker's notation: 

if x ^ y 



y] = 5. 



1 if x = y 



Of course the Iverson's symbol [x — y] is always equal to the Kronecker's delta S X y 
but it has some typographic advantages if x and y are complicated expressions with 
many upper and lower indexes. 

3.4. Random matrices R. 

Definition 2. IfV is a finite dimensional Hilbert space with an orthonormal basis 
ei, ■ ■ ■ , edim v then a hermitian standard random matrix over V is a random vari- 
able M with values in End(V ) such that the joint distribution of the complex matrix 
coefficients My = {ei,Mef) is Gaussian, My = Mji, all My have mean zero and 
the covariance is given by 

(7) E[Mij M u ] = E[MijM^\ = [i = I] [j = k] . 

Alternatively one can define a hermitian standard random matrix (Mij) by saying 
that the following random variables: Ma for all indexes i, ^tM^j, QMij for i < j 
are independent real Gaussian variables with E(Mjj) = for all indexes i,j and 
EM? = ■ s J r7 for all values of index i and E(3?M y -) 2 = E(3My) 2 = 2 di ^ v for 
all i < j. The entries Mij, where i > j are defined by the hermitianity condition 
M (J = Mji. 

One can show that both definitions do not depend on the choice of the orthonor- 
mal basis of V . 

For each A C {1, . . . , N} let us consider a family of hermitian standard random 
matrices i^ Ar ' ,yl ' A1 <= End(7i,4) indexed by /i € 971 such that the entries of different 
matrices are independent. We define a family of random matrices R( N ^ A '^ by 
RW ,A llt = ] a{R [ n )^ i wA ,) g End(ttW) 

where l n A> : TL A — ► TL A denotes the identity operator. 

Intuitively speaking, a matrix R( N )> A >v consists of d N ~\ A \ copies of a <i' A ' x c?l A l 
standard hermitian random matrix. 

As one can see, matrices R A '^ are hermitian and the joint distribution of their 
entries is Gaussian, but different entries need not to be independent. We have: 



R A tfl = r a,^ E[R A ^} = 0, 
and from (|?]) it follows that 



(8) EIR^R^} = E[R A ^R B '^ 



lk 



r r / TT b'r — lr] [jr — k 
OABO^ || 



\r£A 




jr] [kr 
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3.5. Tensors T. The formula (0) can be written shorter if we introduce for all 



AC {I,..., N} and 1 <r < N tensors T^ r kl as follows 



T, 



ij t kl 



ij,kl 

a[i = *][?' = *] 
[i=j][k = l] 



if r e A 
if r A 



We define 
(9) 



-ij,ki 



A,r 

i r j r ,k r l r ' 



what with a small abuse of notation can be written as 



(10) 
Then 

(11) 



rpA _ rpAA 



T 



A.N 



can be written as 



3.6. Examples. First of all note that for the trivial case d — 1 all Hilbert spaces are 
one dimensional and all random matrices R A '^ are in fact scalar random Gaussian 
variables. 

Furthermore, the random matrix ijW^^ is simply a scalar real random Gauss- 
ian variable multiplied by an identity matrix. The random matrix fjWiMil 1 !-^} 
is a hermitian standard random matrix from Definition ^|. 

There is a correspondence between sequences i = («i, . . . , ijv) such that < 
H, ■ ■ ■ , In < d — 1 and the set of integer numbers {0, 1, . . . , d N — 1} given by the 
digit representation of natural numbers in the system with base d: 

i = (ii, . . . , i N ) h-> ii + di 2 H + d N ~ 1 i N . 

. gdN-i of H'^ indexed by 



Therefore we can introduce an orthonormal basis go 
integer numbers: 



.9ii+rfi 2 H l-d N - 1 i N — lii 



f lN ,...,iff) ' 



where < ii, . . . , in < d — 1. In the following, if we want to write an endomor- 
phism M € E,nd(Tt( N )) as a matrix (Mij) 0<i j< (i »_ 1 we shall do it in the basis 

(9i)o<i<d N -l- 

For d = 2 and N = 2 the matrices RW,n,A 

are of the following form: 





O l 








Oio 


an 














°oo 


OOl 








aw 


an 







boi 








&00 







bw 





&n 








bio 





bn 



where 



a oo 
aio 



ctoi 
an 



and 



boi 
bn 



are standard hermitian random matrices. 



Entries of the first matrix are by definition independent of the entries of the second 
matrix. The index fi was omitted, however it should be understood that for different 
values fx the entries of matrices are independent. 
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For d = 2 and N = 3 we have: 



coo 


coi 




















ClO 


Cll 


























coo 


coi 




















ClO 


Cll 


























coo 


coi 




















cio 


Cll 


























coo 


coi 




















cio 


Cll 



doo 





doi 




















c^oo 





doi 














d w 





dn 




















d w 





dn 


























doo 





doi 




















doo 





doi 














dw 





dn 




















dio 





dn 


coo 











eoi 














eoo 











eoi 














eoo 











eoi 














eoo 











eoi 


eio 











en 














eio 











en 














eio 











en 














eio 











en 



where again (c piJ )o<p, 9 <i, (d pq )o< P ,q<i, i e pq)o< P ,q<i are standard hermitian random 
matrices. 

Furthermore 



#{1,2} = 



/oo 


foi 


/02 


/03 














fio 


hi 


/12 


/l3 














/20 


hi 


ill 


/23 














/30 


/31 


fzi 


/33 


























/oo 


/01 


/02 


/03 














/10 


/11 


/12 


/l3 














/20 


/21 


/22 


/23 














/30 


/31 


/32 


/33 



#{1,3} 



ffoo .9oi 
5io 5n 



520 921 
530 331 








9oo .9oi 
5io 3n 



920 921 

930 931 



902 903 
912 913 



922 923 
932 933 



902 903 
912 913 



922 923 
932 933 
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^00 





h i 





h 2 





h 03 













^01 





h 02 





^03 


h w 





hn 





^12 





^-13 








hw 





hn 





^•12 





hi3 


h 20 





/l21 





h 2 2 





h 2 3 













^21 





^22 





h 23 


h-30 





/l31 





/l32 





h33 








fl30 





h 3 i 





h 3 2 





^33 



#{2,3} 



where (f pq )o< P , q <3, (g Pq )o< P , q <3, (h pq )o< P , q <3 are standard hermitian random ma- 
trices. 

3.7. The case of a general covariance r M „. By a small change of definition of 
the matrices R we obtain a more general case. 

Let T^i, be a real positive definite function. For every AC {1, . . . , N} we consider 

a family of random (non hermitian) matrices i?^'" 4 '' 1 £ End(7i y i) such that for 
each pair of indexes i, j we have that the joint distribution of (-Ro ^'Oij ^ £ 9JI is 
Gaussian, 

E(^' A '")y = 0, 

the covariance of real and imaginary parts are defined by the function T: 

1 



Eimi NhA %M4 N) ' A ' u h] = n^(4 N) ' A '%MR ( N) ' A ' u h} 



2d\ A \ 



v 



and the real and imaginary parts are independent: 



E[3?(4 A,) ' A ' M ) ij Q(i? ( 



ij] = o. 



For different choices of sets A or a pair of indexes (i,j) the random variables 
(i?Q W ''" 4 ' A1 )ij should be independent. 
We define hermitian random matrices 



(N),A,fj, 



R 



(N),A,v 




(R 



(N),A,fi 



Note that for the simplest choice of a positive definite function r M „ = 5^ u this 
definition of random matrices jj( Ar ' ,yl,A1 coincides with the definition from Subsect. 



3.4 



Similarly as in Subsect. 3.4 we define 

R W,A,„ = ~ jA{R (N),A,^ 



The joint distribution of entries of hermitian matrices r( ns )' A >^ is Gaussian and 

E[R A <% = 0, 



(12) 



E[R^R^] = E[R^R^} = ^r^T^ kl . 

4. The Main Theorem 
We define a family of random matrices indexed by (i eWl 
S (N)^ = af ] R^ A ^ 

AC{1,...,N} 

where trW is a real- valued function on the set of all subsets of {1, ... , N}. 
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Matrices are hermitian and the joint distribution of their entries is Gauss- 
ian. Alternatively one can define these matrices by giving the mean and the covari- 
ance of the entries: we have 

(13) E[S£] = 0, 

(14) E[S$S& = E[S$SQ = 5> A ) a Tfl )kl . 

A 

In the following theorem we show conditions which the sequence of functions 
(ct^ n >) need to fulfil. Since these conditions may seem quite disgusting, we would 



like to give some hope to the reader by pointing to the equation (15), which gives 
a simple example of a covariance function fulfilling all assumptions. 

Theorem 1. If for each N G N we have that erW is a real-valued function on the 
set of all subsets of {1, ... , N} such that: 

1. (normalisation) for each N G N we have 

E («tf°) a = i> 

AC{1,...,N} 

2. (triple coincidations are rare) 



A 1 ,A 2 ,A 3 Q{1,...,N} 

(distribution of coincidations) there exists a sequence (pi)i>o of nonnegative 
real numbers such that ^2 i>0 p% — 1 and for any k £ N and any nonnegative 
integer numbers riij, l<i<j<kwe have 

N^oo 



e (4 I ) ) 2 -->2 ) ) 2 = n * 



l 1 ,...,A k Q{l,...,N} \<i<j<k 

| A^HAj \ = nij for any l<i<j<k 

4. for each n G N 

jim £ «') 3 -(^ ) ) 2 lPI „^u...u,,,i ^. 

Ai,...,A„C{l,...,AT} 

Then for q — YliLoPiep* we have 



lim E[tr • • • S^'"-] = £ 9^1^ • • • r McmMdm = r[G Ml • • • 



/or every n £ N and fjbx,...,ji n G 971, where G Ml , . . . , G^ ln are q-deformed Gaussian 
variables with covariance T. 

Before we prove this theorem we would like to make some remarks and state 
auxiliary lemmas. 

Remark 1. For a given function crW we define a measure p( N ' on the set of all 
subsets of {1, . . . , N} by assigning to set A the weight (o~^ ) 2 . Then the first three 
assumptions of the theorem can be reformulated as follows: 

1. for each N G N the measure p^ N ^ is probabilistic, 

2. if for each N G N we have that A± , A% , A3 are independent random 
variables with distribution given by the measure p^ N ' then the probability of 
the event: A^ (1 n A3 7^ tends to as N tends to infinity, 
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3. let k be a fixed natural number, if for each N S N we have that , . . . , A^ 
are independent random variables with distribution given by the measure 
then the joint distribution of the random variables \A, H Aj\, 1 < i < j < k 
tends to a product distribution as N tends to infinity. The limit distribution 
of a single variable \Ai n Aj \ is given by 

lim P(\A ( i N) nA i j N) \=k)= Pk . 

Remark 2. The assumption |^ follows from other assumptions, however the proof 
of this is rather technical and we omit it. 

Remark 3. We would like to point out an interesting informal connection between 
stochastical properties of Gi regarded as large matrices and their entries. 

Let us consider Y ^ v = 8^ v . For q = 1 we have that is a family of independent 
Gaussian variables and for q = we have that G^ is a family of free non commu- 
tative random variables. Freeness is an analogue of classical independence; we can 
expect therefore that for in a general case — 1 < q < 1 variables G M are independent 
in some generalised way. 



However, it was proven by Speicher [3p3] that there are only three generalisa 



Hons of the notion of independence of random variables to the non commutative 
setup which would satisfy certain natural properties. These three generalisations 
are: classical independence, free independence (freeness) and boolean independence. 
Therefore except the cases q £ {0, 1} which correspond to the free and the classi- 
cal situation respectively we cannot formally say that the non commutative random 
variables G M are independent in some sense. 

We can of course weaken Speicher 's axioms and treat this "independence" on a 
very informal level. It is worth pointing out that a family of "independent" vari- 
ables Gi is asymptotically represented as random matrices such that the entries of 
different matrices are classically independent random variables. 

Similarly, for the choice of the covariance function T ts = min(t, s) for t, s > we 
obtain a non commutat ive st ochastic process Gt which can be regarded as some kind 



of a Brownian motion | BS1 | and Gt can be asymptotically represented as a matrix 



valued stochastic process. Entries of this matrix are classical Brownian motions. 

Remark 4. The assumptions of the theorem are fulfilled for the following impor- 
tant examples of the functions a: 



Proposition 1. For every real c > the sequence of functions defined by 

/ \ \A\ / x N-\A\ 

fulfils the assumptions of Theorem [J with the sequence p k — yc 2k e c being the 
Poisson distribution with parameter c 2 and 

-(l--4-)c 2 

q — e y d?' . 
In this case the covariance ( fL^ ) takes a beautiful form 
(15) E[SgSEJ = E[S$SQ = 



iv n 



c [ir-lr}\jr-k r } + /_C . r=jr][kr = lr] 



N d V Jn 
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Proof of this proposition will be presented in Sect. |(| 

Proposition 2. For every real number c > the sequence of functions <r ( - Ar - ) defined 
for N sufficiently large by 



((7 W )2 f (^ry if\A\ = [cVN\ 

I otherwise 

where \_x\ denotes the integer part of a real number x, fulfils the assumptions of 
Theorem [J with pk — -gc 2k e~ c and 

Since proof of this proposition is similar to the proof of Proposition |l| we skip it. 

Lemma 1. For any pair partition n — {{ci, di}, . . . , {c m , <i m }} of the set 
{1, . . . , 2m} and any sets A\, . . . , A 2m Q {lj ■ ■ ■ i N} we have 

e n T^ ic „ +1 . d „ id „ +1 < i. 

If furthermore A Ci n A Cj n A Ck — for alll<i<j<k<m then 

E IT T i^"i^+i,i^i<i„+i = Yl ,2|A Ci nA c .| • 

ii...,i 2 ™ Kv<m i<i<i<m " * 1 

lines {c^diMc^dj} cross 

The proof of this lemma will be presented in Sect. ^. 

The following lemma states a well known property of the Gaussian distribution. 

Lemma 2. If the joint distribution of random variables (Xk) is Gaussian and 
E[X k ] = then 

E(X 1 ---X 2m -i)=0, 
E(Xi • • • X 2m ) = J2 nx cl X dl ) ■ ■ ■ E(X Cm X d J, 

where the sum is taken over all pair partitions it = {{ci, di}, . . . , {c m , d m }} of the 
set {1, . . . , 2m}. 

With this preparation we are able to start the proof of the main theorem. 

Proof of Theorem ^|. In the following the sums over ir are taken over all pair parti- 
tions 7r = {{ci, di}, . . . , {c m , d m }} of the set {1, . . . , 2m} and sums over A\, . . . , A n 
are taken over all subsets of the set {1, . . . , N}. Products over v are taken over 
1 < v < m. 

From Lemma ^ follows that for any m £ N and indexes pb%, . . . , \L 2m G 9JI we 
have: 

Etr (SW,Ml...sW,M2m-l) =o 

and furthermore 
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We define 



7T A 1 ,...,A 2 



AN) 



n 



l<i<j<m 
lines {ei.diMcj.dj} 



d 2|A Ci nA Cj .| 



Lemma [j] shows that the corresponding summands in the definitions of U and V 
are equal unless there are some indexes l<p<q<r<m such that A c n A c<! n 
A Cr =/= 0. There are (™) choices of these indexes and again from Lemma [j] and the 
assumption [l] we have 



\uM-vW\<c(fl E E (*2?) a -(*2?) a 

V / A 1 ^ 2 ,A 3 C{1 A 4) ...,A m C{l,..,iV} 



C" 



E (4?)^)Vi?) 2 , 



A 1 ,A 2 ,A 3 C{l,...,JV} 



where C = max r M jM and therefore from the assumption we have 



lim |[/W-yW| = . 

TV— >oo 



We have 



^-sfn^J e 

7T \ D / Ai,...,A 



n 



l<i<j<m 

i> d i}A c n, d n 



1 



>2|A C nA c | 



=Efn r ^)/ n ^7^K i)5 

7T \ U / ^ l<i<3<m 

linos {ej.djMcj,.^} cross 

where measures (A^)) are defined on the set of all sequences (w.y)i<i<j< m , riy S 
{0, 1,2,...} by condition 

a w ({K)})= E (-i?) 2 ---^) 2 . 

A 1 ,...,A fc C{l,...,N} 
jl i nA,-|=ny for all l<i<j<m 

From the assumption H] follows that this sequence converges pointwise to the product 
measure defined on the atoms by 

Ki<j<m 

bmce measures A^) and the measure A are probabilistic, this convergence is 
uniform and the statement of the theorem follows. □ 
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5. The Almost Surely Convergence 

As a simple corollary of Lemma || we have 

Lemma 3. If the joint distribution of random variables X\, . . . 1 X2 m is Gaussian 
and E[X k ] = for all 1 < k < 2m then 

E[Xi ■ ■ ■ X 2m ] - E[Xi ■ ■ ■ X m ] E[X m+ i ■ ■ ■ X 2m ] = 

=^2 n E i x cM, 

7r' l<v<m 

where the sum is taken over pair partitions tt' — {{ci, d±}, . . . , {c m , d m }} of the 
set {1, . . . , 2m} which have additional property that there exist i£ {1, . . . , m} and 
y G {to + 1, . . . , 2m} such that {x, y} G w' . 

We define a permutation a : {1, . . . , 2to} — * {1, . . . , 2to} by er(fc) = k + 1 for 
k (ji {to, 2to}, ct(to) = 1, <t(2to) = to + 1. 

Lemma 4. Let tt' = {{ci, di}, . . . , {c m , d m }} be a pair partition as in Lemma 
i.e. for some 1 < k < to we have Ck G {1, . . . , m\, dk G {to + 1, . . . , 2m}. 
Then 

— V TT t Acv 1 

MN Z—l 11 i^i^^i'M"^ 



d 2N 

i 1 i 2m 1 < t> < 7ri 



This lemma follows directly from Lemmas ^| and |t] from Sect. |(]. 



< 



Proposition 3. 

Var[tr S {N) ^ ■ ■ ■ S {N) ^} < 

<r V (.t (n) ) 2 <.t {n) ) 2 1 

- 2^ y (T A 1 ) { - (T A m ) d 2\ Al \(A 2 U---UA m )\' 

A 1 ,...,A m C{l,...,N} 

where C = (2m)!! max |r Mp , M J . 

Proof. We define /J. m +k = Mfc- I n the following sums over 7r are taken over all 
pair partitions n = {{ci, d\}, . . . , {c m , d m }} of the set {1, . . . , 2m} with additional 
property that there exist x G {1, . . . , m} and y G {m + 1, . . . , 2to} such that {x, y} G 
tt'. 

From Lemma we have that 

Var[tr S^'" 1 • • • S^*-] = 
= E[(tr SW-" 1 • • • SW^) J ] - [E(tr S^-" 1 • ■ ■ S (JV ^'")] 2 = 

~~ TfilV I -Ul'L»J. li2 • 0. m _ij m 0. m ji 0. m+1 . m+2 • <3j2m-lJ2m'Jj2mJm+lJ 

° il,...,i 2 ™ V 

-E[S& • • • Sfcr^Sfo] E[S£ +1 . m+2 • • • S&Z\ i2m S& im+ i}) = 

tt' i 1 i 2 ™- Kv<m 



T W JN) 

d 2N 

it' Ai,...,A 2m C{l,...,iV}ii,...,i 

X IT ^ A o nj A d ^fi Cv ,fJ, dv T i a v - ia (c v ) t ld v la(d v ) < 

1< v < m 
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where in the last inequality we used Lemmas |^ and ^. □ 

If {Mij)\<i.j<K is a hermitian matrix with eigenvalues Ai, . . . , Xk we define a 
probabilistic measure vyi on the real line R by the following 

l<n<K 

Theorem 2. // the assumptions of Theorem are fulfilled then there exists an 
increasing sequence of natural numbers (Ni) such that the sequence of measures 
Vg(Ni) almost surely converges weakly to the measure v q . 

Since for < q < 1 the support of the limit measure v q is compactly supported, 
the theorem follows from the following stronger statement. 

Theorem 3. // the assumptions of Theorem [I] are fulfilled then there exists a se- 
quence (Ni) such that the limit 

lim tr ^'1* • • • S im) '^ = r(Gu 1 • • • G M J 

rn — >oo 

holds almost surely. 

Proof. Our goal is to construct a sequence (Ni) such that 



5]E [(tr S W,w... S (^ m _ T ( G(ii ... Gji j) ; 



< oo 



holds. 

However, 

E 



(tr S<- N ^^ ■ ■ • S^'-"™ - t(G Mi ■ ■ ■ G^JY 



(E [tr S^'" 1 • • • S (N) ^ m - r(G Ml • • • G„J]) + Var [tr S^'" 1 • • • S^'""*]. 

The first summand converges to by Theorem [l]. 
From Proposition || we have that 

Var[tr S^'" 1 • • • S^'"™] < 

<r V (.-t [n) ) 2 (.t [n) ) 2 1 

- 2^ ^ Ai J ^°A m J rf 2|Ai\(A 2 U-UA m )| 

A 1 ,...,A m C{l,...,N} 

From assumptions |] and || it follows that this expression converges to 0. □ 

6. Proofs of Technical Lemmas 
Lemma ^ follows directly from the following two lemmas. 

Lemma 5. For every n,M £ N if for all 1 < v < M we have A v C {1, . . . , N}, 

1 < Pv, qv,r v , s v < n then 

(i6) E II T^ iW = II X II t--^,. 

i 1 ,...,!" 1<«<M l<r<N 0<i 1 ,...,j'"<d-l l<-u<M 

Proof. This lemma is a direct consequence of Eq. (0) and din) . □ 
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For sets A\,...,A2 m C {1,...,7V} and a pair partition tt = 
{{ci, di}, . . . , {c m , d m }} of the set {1, . . . , 2m} we define 

(17) & Ai,...,A 2m ,n = I Y Y\T Ac -' r ■ - 

0<j 1 ,...j 3m <<J-l u 

1 J = j^+i]^ = jdu+i] if r £ Acu 

d 2-, 11 | lyc. •«w][ f .+i =i *. ] ifreA Cu ■ 

0<J j-")J 2m <d— 1 l<«<m v " 

It should be understood that j 2m + 1 = j 1 . 

Lemma 6. For any Ax,...,A2 m C {1, ...,./V} cmc? a pair partition n — 
{{ci, di}, . . . , {c m , d m }} we have 

(is) o< [| e^i.-^am."- < 1. 

l<r<W 

If furthermore A Ci n A Cj n A C(c = /or alll<i<j<k<n then 

d9) n e^-.^-= n -^W 

Kr<N l<i<j< m a * 3 

linos {ci.djMCj.dj} cross 

Proof. Let us consider an expression of the following type 

(20) £ [f 1 = Z 1 ]-- Lf 2m =i /2m ]- 

0<3 1 ,...,i 2m <d 

We can represent this expression by a graph (see for example ]Do| ) with 2m vertices 
which are labelled by variables j 1 , . . . , j 2m and with vertices j Ei and j^* connected 
by an edge for all 1 < i < 2m. We see that the nonzero summands in (|2^) come from 
indexes (j 1 , . . . ,j 2m ) such that to all vertices in the same connected component of 
the graph is assigned the same value. Therefore the expression fl20| ) is equal to d M , 
where M denotes the number of connected components of the graph. 

Let us fix an index r. We shall apply the above observation to evaluate 
QAi,...,A2m,ir 4 L e t n be the number of indexes v such that r € A Cv and let v%, . . . , v n 
be all such indexes. We consider the following cases. 

1. If n — then we obtain a graph of type presented in Fig. [j]. This graph has 
one connected component and therefore 0^i>—>-A2m,w — \ 

2. If n = 1 then we obtain a graph of type presented in Fig. ^|. This graph has 
two connected components and therefore ® Ax ' ~ > A2m >' K = 1. 

3. Suppose that n = 2. If the lines {c Vl ,d Vl } and {c V2 ,d V2 } do not cross (Fig. |^ 
and ^) then corresponding graphs have three components and and therefore 
Q A u ...,A 2m ,* = i. 

But if the lines {c Vl ,d Vl } and {c V21 d V2 } cross (Fig. ||) then the graph has 
only one connected component and therefore Q A i>- = -L . 

4. In the general case n > 3 the graph has at most n + 1 components and there 
is a factor therefore < 9j5^-. A 2m,» < i 

Therefore if for a given partition tt sets A Cl , . . . , A Ck have a property that the 
common part of each three of them is empty then 



(2i) n e^--^= n 



i 



l<r<7V l<i<j<k 

linos {o i ,d i },{o j .d j } 



d 2\A Ci nA c .\ ■ 




Figure 1. 




Figure 4. 



In the general situation the expression Yl 1<r<N Qr 1 ' "^ 2 ™'* * s a rea l number 
from the interval [0,1]. □ 



2(1 



PIOTR SNIADY 




Figure 5. 



We recall that the permutation a : {1, . . . , 2m} — > {1, . . . , 2m} was defined by 
er(fc) = k + 1 for fc ^ {m, 2m}, cr(m) = 1, a(2m) = m + 1. 

For sets A\,...,A2 m Q { 1 , . . . , A^} and a pair partition tt = 
{{ci, di}, . . . , {c m , d m }} of the set {!,..., 2m} we define 



(22) 



c u cr(c u ),d w cr(d tJ ) 



0<j 1 ,...,j 2m <d— 1 l<ti<r 



e n 

0<j 1 ,...,j 2m <d-l l<v<r 



■■ j dv ] if r S ^4 Cl , 



ir„-c„ 



<x(d„)j^ CT (c„) 



Lemma 7. For any Ai,...,A2 m C {1, ...,iV} arte? a pair partition n = 
{{cx,di},...,{c m ,d m }} we have 

(23) < T^ 1 '- ^ 2 ™^ < 1. 

If furthermore r G \ Uj^fc A? ^ e7t 



(24) 



•>Ax,...,A2 m ,ir _ 



d 2 



Proof of this lemma is very similar to the proof of Lemma |^ and we will skip it. 

Proof of Prop. [J. The measure p^ N ^ on the set of all subsets of {1, . . . , TV} such 
that for any AC{1,..., N} we have 



pW({A})=(S=) U-^= 



N-\A\ 

'n ) V Vn, 

is simply a Bernoulli distribution and the assumption [l] of Theorem |l| is obviously 
fulfilled. 

Therefore if A\, A%, A% are independent random variables with distribution given 
by then 

P{lo : Ai(w) n A 2 {lo) n A 3 (lu) ^ 0} < 

. 3 

< P{lo : k G Ai(w) n A 2 (w) n A 3 (w)} = N 

l<k<N 

and the assumption || of Theorem |l] is fulfilled. 

Let Ax, ... , Ah be independent random variables with distribution given by p^ N ^ 
and let us consider the probability of following event: \A p D A q \ — n pq . The 
discussion from the previous paragraph shows that this probability is equal (up to 
an error of order 0(iV~ 5 )) to the probability of the event: \A P H A q \ = n pq and 
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furthermore A p n A q n A r = for all p < q < r. We shall now evaluate the latter 
probability. 



There are , j r - — n/^'_y - — yj choices of disjoint sets (i?pg)i<p< q <fc such that 



|-Spg| — n pq- 



For i € B pq the probability of the event: i G A p D A q and £ ^ A r for r ^ {p, q} 
is equal to 

For £ ^ B pq the probability of the event: £ £ A p n A g for all 1 < p < q < k 
is equal to 

)k / \ fe— 1 



Therefore the probability of the considered event is equal to 
N\ 



1 



)fc / \ fe— 1 

, c / c 
+ fc^ l--7= 



where all sums and products are taken over 1 < p < q < k. 

After short calculations one can show that the limit of this expression as N tends 
to infinity is equal to 

1 



n 



n pq \ 



c 2n pi e 



l<p<q<k 

and therefore the assumption || is fulfilled. 

The assumption [| follows from the observation that the distribution of the ran- 
dom variable \A± \ (A% U • • • A n )\ is binomial and simple computations. □ 
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